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Abstract. We discuss freeze out on the hypersurface with time-like normal vector, trying to answer how 
realistic is to assume thermal post freeze out distributions for measured hadrons. Using simple kinetic 
models for gradual freeze out we are able to generate thermal post FO distribution, but only in highly 
simplified situation. In a more advanced model, taking into account rescattering and re-thermalization, 
the post FO distribution gets more complicated. The resulting particle distributions are in qualitative 
agreement with the experimentally measured pion spectra. Our study also shows that the obtained post 
FO distribution functions, although analytically very different from the Jiittner distribution, do look pretty 
much like thermal distributions in some range of parameters. 
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1 Introduction 

Since the very beginning of heavy ion collision physics 
fluid dynamical models were widely used for simulation of 
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such reactions. Their advantage is that one can vary flexi- is an ideahzation of a layer of finite thickness (of the or- 

bly the Equation of State (EoS) of the matter and test its der of at least mean free path or collision time) where the 

consequences on the reaction dynamics and the outcome, frozen-out particles are formed and the interactions in the 

This makes fluid dynamical models a very powerful tool matter become negligible. The dynamics of this layer is 

to study possible phase transitions in heavy ion collisions - described in different kinetic models such as Monte Carlo 

such as the liquid-gas or the Quark-Gluon Plasma (QGP) models pillUj or four-volume emission models jlllll2lll3l 

phase transition. For highest energies, achieved nowadays I14II15| . Under such an assumption FO can be pictured as a 

at RHIC, hydrodynamic calculations give a good descrip- discontinuity in rclativistic flow, where the kinetic proper- 

tion of the observed radial and elliptic flows ^|2]|3lE]|ol ties of the matter, such as energy density or flow velocity 

E], in contrast to microscopic models, like HIJING [7| and change suddenly. 
UrQMD [HI. 

It took a long time to develop general theory of discon- 
Freeze out (FO) process is an important and neces- tinuities in the rclativistic flow. The story goes like this. In 
sary part of the fluid dynamical modeling. Generally say- 1945 Taub [IBj discussed discontinuities across propagat- 
ing, fluid dynamical models are determined by the initial ing hypersurfaces with space-like normal vectors. If one 
state, EoS and the way FO is realized, in addition to well applies Taub's formalism to FO hypersurface with time- 
known fluid dynamical equations. FO describes the last like normal vector, one gets a usual Taub adiabat, but the 
stages of the reaction, when the matter becomes so di- equation of the Rayleigh line yields imaginary values for 
lute and cold that particles do not interact anymore, their the particle current across the front. Thus, these hyper- 
momentum distribution freezes out, and particles stream surfaces were thought unphysical. However more recently, 
freely toward the detectors. FO models allow the evalu- Taub's approach has been generalized to these hypersur- 
ation of two-particle correlation data, different flow an- faces as well [T7| and the imaginary particle currents aris- 
alyzes, transverse momentum and transverse mass spec- ing from the equation of the Rayleigh line were eliminated, 
tra and other observables. There are many different theo- Thus, it is possible to take into account conservation laws 
retical approaches for the freeze out problem. Sometimes exactly across any surface of discontinuity in rclativistic 
these approaches are oversimplified to the extent that con- flow. The corresponding equations are 
servation of energy and momentum is violated (see Ref. 

miEa for an overview). T^:fJ<J, = T^^d^. , (1) 

A basic standard assumption in fluid dynamical mod- 
eling of heavy ion collisions is that freeze out happens 

across a hypersurface in space-time. Such a hypersurface UpreUpre^'^tJ. — ''^postUpostd'^ tJ- i (2) 
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where da^ is the unit 4- vector normal to the discontinuity aUstic four- volume emission model then final distri- 
hypersurface, butions have more complicated form, although they are 



Tf"" ^{e+p)ut'u'' ^pg'"' (3) 



cut-off as required. Unfortunately, so far we do not have 
good analytical expression to fit post FO distributions for 
is the energy-momentum tensor, the hypersurface with space-like (or, generally, any type 



Nt" = nu'' (4) 



of) normal vector. Therefore in most calculations scien- 
tists prefer to model FO based on time-like hypersurfaces. 
is a baryon current (usually in heavy ion collision mod- For the ultra-relativistic heavy ion collisions this is justi- 



eling we neglect electric charge current) These consist of hed by the Bjorken model, where the natural choice of FO 

local thermodynamical fluid quantities: the energy den- hypersurface is t = tfo- 

sity e, pressure p, baryonic density n and the collective ^his paper we ask the following question. Well, we 
four-velocity u'^ = -^^=^{l^v). know now that FO through the hypersurface with space- 
Since 1974 people were using the Cooper-Frye formula jj^e normal vector leads to non-equihbrated post FO dis- 
UHl to calculate final particle spectra. Much later it was re- tributions, but can we really assume thermal post FO dis- 
alized that the Cooper-Frye FO description has a concep- tributions if the FO hypersurface is time-like? This is not 
tual problem of negative contributions, when it is applied restricted by general theory, but how realistic such an as- 
to the FO hypersurface with space-like normal vector [H gumption is? From the experimental data we know that 
PD| | 2H22 ||23|| 24 | |^Sl Eni -^ This leads to the non-equilibrated pj^j^ g^j^j charged hadron (which are actually pion domi- 
cut-off post FO distributions: nated) transverse mass spectra both at SPS and at RHIC 

|29ll3l)j are not the Jiittner type - the slope parameter is 



decreasing with increasing pt. 



fFo{r,p;T,n,u-',da'^) ^ fFoir,p;T,n,unO{p^ da,,) . 

(5) 

The simplest example of such cut-off distribution is the We generalize the simple model presented in Refs. ^ 



cut Jiittner distribution 28 , proposed in Ref. In EHHl for the FO hypersurface with timc-Uke normal 

Refs. pni^^lEHll^ it was shown that one can obtain cut sector (the model is also not ID anymore). The philosophy 

Jiittner post FO distribution in a oversimplified kinetic oi this work is also similar - first we will show that in the 

model. If one improves this model, taking into account oversimplified model the obtained post FO spectra are 

rescattering and re-thermahzation, or if one uses more re- ^ally thermal ones. Then in more reahstic model, taking 



^ To our knowledge the first attempt to modify post FO dis- into account rescattering in the still interacting gas, we 
tribution has been done by Yu. Sinyukov in 1989 |57|. For the will obtain post FO distribution in the more complicated 
critical review of this work see [53 . non-equilibrated form, which is in qualitative agreement 
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with experimental data, mentioned above. Actually the (below we will use shorter notations, //reelT) and /^^(t), 
main goal of this work is not quantitative calculations, but but the meaning is the same). At the initial moment, t = 
a qualitative understanding of the phenomenon in simple tq, the distribution ffree vanishes exactly and /i„t(ro) = 
models. fjuttneriTo,r,p,T,n,u^). In this work we will use an ex- 

tended definition of the classical Jiittncr distribution , 
2 Simple freeze out model which is adequate at ultra-relativistic energies and allow 



As it was mentioned above our basic model is in some 
sense similar to the one discussed in Refs. 
but on the other hand it is much more realistic and also 
habitual for the people working in ultra-relativistic heavy 
ion collision field. Instead of very special geometry (ID, 
infinitely long tube with particle source on the left end 



for particle creation. Thus the particle density, which de- 
pends on T, is not the same as the conserved baryon 
charge density, n. Then, fi„t gradually disappears and 
ffree gradually builds up as t increases. The most simple 
kinetic model describing the evolution of such a system is 
the following: 



and vacuum on the right one) used in Refs. |21II22II23II24| . drfint{T) = — /mi(''') , 

this model can be combined with any relativistic hydrody- 

namical model. Below we will use proper time coordinate (6) 

keeping in mind the Bjorken hydrodynamical model, as ^ 

drffree{r) = — /i„t(r) , 

the simplest but still reliable one, although the similar ''7 

treatment could be done for any pre FO state. where Tf is a parameter with dimension of time deter- 

Let assume that at some moment t — tq our matter mining how fast particles freeze out. It is of the order of 

is still completely equilibrated and described by the ther- collision time. In general case Tf may depend on particle 

mal distribution fjuttneriTo,r,p,T,n,u^). Starting from density and momentum of the particle, but in our simple 

this moment matter starts to gradually freeze out. At model we neglect such effects and consider Tf — const. 

each particular moment the FO hypersurface will be Our assumption becomes particularly bad for the large 

T = Ti, and we will follow this hypersurface, i.e. we will times, when the particle density is very low - we will dis- 

describe our system from the reference frame of the front cuss this problem separately in section|^31 If this approx- 

(RFF), where normal vector is dcr^ = (1,0,0,0). We can imation is not done, i.e. the relaxation time is assumed 

describe the FO kinetics of such a system assuming that to be momentum and density dependent, that would be 

we have two components of our distribution (12II13II14L a large step towards the full solution of the BTE, which 

I21II22II23II24) : ffree{T, r,p), describing frozen out particles, is then a much more involved problem, and in most cases 

and /i„t(r, r,p), describing particles, which still interact cannot be performed analytically. Nevertheless even this 
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strongly simplifying assumption leads to nontrivial results term and the re-equilibration. If we include the collision 

as we shall see later. terms explicitly into the transport equations ((HJ this in 

The interacting component, fint, will gradually freeze general case leads to a combined set of integro-differential 

out in an exponential rarefaction with time. equations. We can, however, take advantage of the relax- 



fintiT) — fjuttneriTa)e . (7) 



ation time approximation to simplify the description of the 
dynamics. In this framework the two components of the 
Then momentum distribution develop according to the modified 



ffreeir) = fjuttneriro) [l - e . (8) 



coupled differential equations 



,2 



A. J- . -u J- -,1. J. •• 1 drfuAr) ^ -^fintir) + [feq{T) - fintir)]^ , 

At T — > oo the distribution J free will tend to the original ~f ~r 

Jiittner distribution, while fi„t will vanish. The frozen out ^g^ 

particles keep their velocities and consequently the volume ^ 

drffreeir) = — /i„t(r) . 

occupied by the frozen out matter will gradually increase 

with time Thus, the interacting component of the momentum 

This is a highly unrealistic model, similar to the one distribution, described by eq. ©, shows the tendency to 

resulting in the cut Jiittner distribution in Refs. fUmmi approach an equihbrated distribution with a relaxation 

121 for the FO hypersurface with space-like normal vector. r,. Of course due to the energy, momentum and con- 

We can improve this model taking into account rescatter- served particle drain, this distribution, fegir) is not the 

ing in the interacting component in the same way as it initial Juttner distribution, but its parame- 

was done in the References above. ^^rs, neg(r), Teqir) and u^,(r), change as required by the 

conservation laws. 



If we do not have collision or relaxation terms in our 
transport equation then the conservation laws are triv- 
ially satisfied, because the change of the full distribution 

function is drf = drfintir) + drffreeir) = 0. If, however, 

distribution with decreased amplitude is, of course, highly . i ^ ^ .^i ^ --u j. 

° collision or relaxation terms are present, these contribute 



2.1 Freeze out model with rescattering 

The assumption that the interacting part of the distri- 
bution, fint, remains after some drain the same Jiittner 



unreaHstic. As suggested in Refs. pTll22ll2:^ll^ rescatter- 



to the change of T'^" and N'^, and this should be consid- 



ing within this component will lead to re-thermalization 2 ^^^^^^ ^hat in the ultra-relativistic case the total number 

and of particles is much bigger than the baryon charge and depends 

re-equilibration of this component. Thus, the evolution on T. Thus the collision term in system @ gives a nontrivial 

of the component f^t is determined both by the drain and nonvanishing contribution. 
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ered in the modified distribution function /i„i(T). fn tliis tensor of this component have the usual forms, eqs. (|3I4|) . 
case from the conservation laws we get: Thus, we can write eqs. ()12I13|I in the following forms: 

dNt„, = -dN^^^^ = [ ^p'^fMr) = nir)u^ir) = n(ro)u^(ro)e-^ (15) 



(10) 



and 



Tf J po Tf 



and 

(e(T)+p(T))u^(T)u-(r)-p(T)5^'^ 



Tf J po""'^' Tf »*' =MTo)+p{To))u^{To)u-'{To)-piTo)gne .(16) 

(11) 

Using the expressions derived in Ref. [21| we can easily 



where we used kinetic definitions of TV and T'"'. The 
above equations can be easily solved: 



and 



find the solutions for ^(t), e(T), p{t) and u^^{t). First of 

all, it is nice to sec that now we do not have a problem with 

^inti'^) — (12) the flow velocity definition (for the FO model discussed in 

Refs. |21II22II23II24 | the Landau and Eckart definitions of 

^-^0 the flow velocity were giving different results), because 

TZ{r)=T^''{ro)e- ^/ . (13) 

the velocity change vanishes: 
2.2 Immediate re-thermalization limit du'i , At) = Z\^''(t) ^^"'^^^ = 

Eckart\ I \ J TI^t) 

Z\^'^(r) dT (t) u'^(t) 

As a flrst approximation to this solution let us assume — du'^^^^^^{T) = r \ n/\ ^ = 0, (17) 



that Tr <C Tf, i.e. re-thcrmalization is much faster than 
particles freezing out, or much faster than parameters, 
neq{T), Teq{T) and u^q{T) chauge. This leads to 

fintir) « feq{T), for < T/ . (14) 

This assumption may be unrealistic for large times, when 
particle density becomes very low. We shall discuss the 



e(r) + P(r) 

where A^'^{t) = g^" — u^{t)u'^{t) is a projector to the 
plane orthogonal to u^{t). Thus, 

u^'(r) = u^(ro) = = const. (18) 

Then, 

de(r) = dTZir) ^ e(r) = e(ro)e . (19) 



problem of the large times in the next section. „ n mh I'. v i . • ; i • r . i i 

° Irom eqs. (115118(1 we obtam the expression tor the baryon 

For feq{T) we can assume the usual Jiittner form at any (;[gjjgj^y 



r with parameters ?t.(t), T(r) and u^(t), where u^{t) is 
the flow velocity in the RFF frame. Under such an assump- 
tion the interacting component is all the time in equilib- 
rium, and therefore baryon current and energy momentum p(t)— p(ro)e ^-f . (21) 



n{T) = n(To)e . (20) 
Similarly, from eqs. P6I18I19|I we obtain 
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This last expression means that from the moment our Now, based on the second equation in system 0, we 
gradual FO has started the pressure in the interacting can find the distribution function of the frozen out parti- 
component is fixed by the model and conservation laws, cles. 



(2^) 



and one does not need in addition the EoS as usual. The 4 

ffreeir) 

EoS during FO looks like: 

p{r) = ^e{r). (22) where 



El ( _4!^e— -El 



Tiro) J V T{to 



(29) 



2.3 Freeze out of massless pion gas 



Ei{y)= / -dy, (30) 

-oo V 



is the Exponential integral function. 



To proceed further we need to make some assumptions 4 / jj/j^ \ 

ffree{T)\r^oo ^ "TTT^i^^ ^ , (31) 



about the type of matter we are dealing with. First of all ^^""-^ ^ T{tq) J 

we assume that our system is homogeneous (this assump- figure [H shows the evolution of the /,„t (upper sub- 

tion will be used for the next sections as well). Let us take P^^*) and f^ree (lower subplot) with time, Ar (r - 

the simplest case of massless pion gas. Then, EoS is ^o)/r/. The initial flow velocity (which is actually kept 



e = osbT^ , (23) 



constant all the time) is w = (inspired by Bjorken model). 
T{tq) — 170 MeV - we want to start our simulation right 
^ ~ 3^ ' ^^^■^ after hadronization. We present the distribution functions 

and it will stay like this during all the FO process. Using vs. absolute value of the momentum, as for zero flow ve- 
eq. (|19|) we find temperature as a function of r: locity all the directions are equivalent. 



T(T)=r(ro)e . (25) 



We would like to underline the following points about 

these graphs. One may notice that for the small times t 

Thus, we can write the expression for the Maxwell-Boltzmann 

(At < 3) ffree looks pretty similar to some thermal distri- 

(MB) distribution function of the interacting component 

bution - a straight line in the logarithmic scale in our case 

(we are working in the system, where h = c = 1 and as- 

- although the analytical expression is very different. For 

sume for simplicity only one type of pions. i.e. degeneracy 

larger Ar values ffree starts to show clear deviations from 
the thermal-like distribution. Its non-thermal behaviour, 



fint(T) = e~ "^t^) = e ^(^o)" . (26) namely the decrease of the slope with increasing momen- 

(27r)'^ (27r)^ 

turn, is in qualitative agreement with the experimentally 

/««t(T)|r^oo ^0, (27) 

measured pion and charged hadrons spectra |29II30| . 

as it should. For such a distribution 

For the very large times our distribution becomes very 

3 

^SB = ^ ■ (28) strongly peaked at low momenta. This is not a surprise, 
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since our temperature is going down exponentially and for where we have used eq. (|28|) . The entropy density for t = 
large times we keep only particles with very low momenta oo, when all the particles are frozen out, can be calculated 
in our interacting component, which keeps gradually freez- through eq. H33|) using the distribution function given by 
ing out and cooling down. expression We have done this integration numerically 

As it was mentioned before our model becomes unre- and obtained for the ratio of entropy densities: 



i?= ^^^—^ = 0.9449, (35) 
s(to) 



alistic for large times, but there is staightforward way to 
solve this problem. After some time the particle density 

will be so low that the particles will have no chance for fur- for any initial temperature, T(to)^. 

ther interaction. So, we have to stop at some moment and Thus, we see that the entropy density is slightly de- 
transfer the rest of still "interacting" particles (probably a creasing during our FO process. Since we have a gradual 
few percents) to frozen out component. In any case, since FO this means, that in order to avoid a decrease of total 
our FO is a very fast exponential process, the uncertainty entropy, our FO model has to be accompanied by an ex- 
generated at large times will be rather small. pansion, which is natural, as mentioned before, and which 

should lead to a volume increase by at least 6% during 

2.4 Entropy condition FO. This does not sound unrealistic. 



To obtain a physically realizable result, we have to check 
the condition for entropy increase: 



2.5 Freeze out of nucleons 



Now we would like to study the situation for baryon rich 
matter - in this case the distribution functions will be dif- 
ferent, since we have to take care about baryon density 
conservation in addition. If our system contains some nu- 
J P j-^^-^ |(27r)'^ /("T")! ^ 1 cleons in addition to pions (for simplicity we assume that 



Sr~Sro>0, (32) 

where, for both the equilibrium and the non-equilibrium 
distributions the entropy current is given as 



(33) 



there are no other mesons in the system), then the situ- 



This condition is not necessary to obtain a solution of the ^tion is the following. The EoS at the pre FO state was 

freeze out problem, but it should always be checked to complicated then eq. (jMl, but for the FO times it is 

exclude non- physical solutions fixed by the eq. it^ . We assume that the temperature is 

In our case the original entropy density is given by the ^'^^ ^nd the chemical potential is not too large - such that 



well known expression for the Maxwell-Boltzmann gas: 



It is not a surprise that R is independent on T{to), since 
for the massless gas we have only one dimensional parameter 



4 e(To) 4 3 

s('^o) = '^T{to) ^ 7r2^ '■^"■^ ' ^^'^'^ ~ ^("^o) ^^"^ correspondingly s oc T^(to) 
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our system is pion dominated (this is reasonable assump- 
tion for ultra- relativistic heavy ion collisions). For such a 
system the temperature will be still basically defined by 
pions and we can connect e and T through equation 123|l 
and correspondingly eq. (|25() is also valid. ^ 

The expression for the distribution function of the in- 
teracting component now becomes (see for example '31'): 

1 



fintir) 



n{T)e "'^(-'^ , (39) 



If we would use Bose and Fermi statistics we could formu- 
late this approach in the following way: 



and 



(36) 



(37) 



From the last equation using eqs. 12512011 . we can find that 

m(^) = fi{ro)e~^ , (38) 

and, thus, fi is decreasing faster than T, and therefore if the 
inequality /i <g T, holds for r = ro, it will be even stronger for 
higher r. 



Thus we obtain 

f^nt{T) - CT-^^\T)n(T)e-'-^ 



= C T-^^^{To)n{To)e ' e ^(-o)- ^ , (40) 

where C = j^;;^^- Again we see that f^nt{T)\r^oo ^ . 

Now, based on the second equation in system (jSJ, we 
can calculate the distribution function of the frozen out 
particles. 



2T:m?TK2{m/T) 

where m is the mass of the nucleon, K2{x) is the Bessel 
function of the second kind. Of course, for the proper 
treatment of the nucleons we would have to use Fermi 
statistics, but in this work we restrict ourselves to MB where 
distributions first of all to get an analytical/semianalytical 
results, and secondly because MB distributions are never- 
theless frequently used to fit experimental data, see for ex- 
ample Pn]. Since for the FO in heavy ion collisions T <C m 
we can use the expansion of the Bessel function: 



(41) 



r{a,x) = I g-^x^-^dx 

J X 

is the incomplete gamma function.^ 

//ree('^)|r— >oo ^ 

4CT-3/2(ro),i(To) {p'-u^ - mf'^ r {~^-,p^^u^ ~ 



(42) 



(43) 

The nucleon distribution functions are presented in 
Figure|21 We choose the initial baryon density much smaller 
than normal nuclear density, n(ro) = 0.0025 jmT^ ^ to jus- 
tify our approach |3J - 

3/2 



= 31 MeV < T(to) = 170 MeV . 

Again ffree looks similar to some thermal distribution. In 

this case f free is saturating very fast - already the ffree 
^ Please note that the Exponential integral function, eq. 

13011 ■ is actually a special case of incomplete gamma function: 

r{Q,x) = -Ei{~x). 
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(At = 3) is a good approximation of the final //ree(Ar = 4 Conclusions 
oo). This is due to the fact that fint is freezing out much 
faster for baryons since, in addition to the exponentiaUy 
decreasing temperature, it now has the baryon density, 
n(T), in front, which drops down even faster than T. 



3 Freeze Out of dynamically expanding 
system 



Our analysis shows that, although the thermal post FO 
distributions are not restricted by general theory for the 
FO hypersurface with time- like normal vector, it may be 
too optimistic to fit particle spectra with thermal distri- 
butions. In our simple kinetic models we are able to cal- 
culate post FO distributions analytically, and we may in- 
deed reproduce thermal post FO distribution, but only in 
highly simplified situation. If we take into account rescat- 
tering and consequently re-thermalization in the interact- 
ing component, the post FO distribution gets a more com- 
plicated, non-thermal form. The resulting particle distri- 
butions are in qualitative agreement with the experimen- 
tally measured pions spectra. 



As it was shown in section IT^ we have to apply our FO 
model to the expanding system, what is natural in heavy 
ion collisions. The simplest, but still nontrivial, choice we 
kept in mind from the very beginning is the Bjorken hy- 
drodynamical model. 

Combination of this FO scenario with Bjorken scenario 

-L • J T f 1 i J- j-1 To clarify the principal place of our work in the field 

IS straightiorward. Let us, tor example, concentrate on the t- t- t- 

J rpi • •;• 1 ; ; • • . wc shall makc a comparisou wlth thc cut Jiittucr approacli 

energy density, ihc initial state is given at some Ti„ - ( = 

I \ T? ,-\,- i J -i u 1 1 i J of Bugaev [191, eq. ISl which is the only other post FO 

e(Tm). l^rom this moment energy density can be caicuiated ° ^ j 

■ ,1 T-,- 1 1 J J -1 1 1 r ii i- distribution published up to recently, which satisfies the 

m the Bjorken hydrodynamical model irom the equation: ^ '- 

FO requirements for space-like FO surfaces. Afterwards, 

de _ e+p 

dr ~ T Bugaev's approach was improved for space-like surfaces 

in Refs. |21ll22j and in subsequent publications by this 
research group. The present work completes this improve- 
ment work and extends it to time-like FO hypersurfaces 
also. Thus, more precisely, this approach is an improve- 
Eq. is valid until moment r = tq, when FO starts, ment compared to Ref. ^\ just as the works [?TII^ . 



For the ideal ultra-relativistic gas with EoS p ~ e/3 we 
obtain: 

>4/3 



s(r)=e(r,„)(^) ■ (45) 



Then during FO process energy density is given by the This work actually completes the approach introduced in 
expression: |21II22| by extending them to time-like freeze-out surfaces. 



4/3 The approach is not the best of all possible approaches, 

e(r) = e(ro)e"— = e(r„0 ( — ) e"^ . (46) 

''o / it includes a few simplifying assumptions for the more 
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transparent presentation, but it is an important advance 9. L.V. Bravina, I.N. Mishustin, N.S. Amelin, J. P. Bondorf 

compared to the cut Jiittner approach. and L.P. Csernai, Phys. Lett. B354, (1995) 196. 

Another important point one can learn from our study 10. L.V. Bravina, I.N. Mishustin, J. P. Bondorf and L.P. Cser- 
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Fig. 1. Figure shows the evolution of the fint (upper subplot) and ffree (lower subplot) with time, At = (t — ro)/r/. The 
initial flow velocity (which is actually kept all the time) is = 0, T(to) = 170 MeV. We present distribution functions vs. 
absolute value of the momentum, since for zero flow velocity all the directions are equivalent. 
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